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The chemical potential of a hard-sphere fluid can be expressed in terms of the contact value of 
the radial distribution function of a solute particle with a diameter varying from zero to that of 
the solvent particles. Exploiting the explicit knowledge of such a contact value within the Percus- 
Yevick (PY) theory, and using standard thermodynamic relations, a hitherto unknown PY equation 
of state, pi pksT — — (9/n) ln(l — rf) — (16 — 31?d /2(1 — n) 2 , is unveiled. This equation of state turns 
out to be better than the one obtained from the conventional virial route. Interpolations between 
the chemical-potential and compressibility routes are shown to be more accurate than the widely 
used Carnahan-Starling equation of state. The extension to polydisperse hard-sphere systems is 
also presented. 
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Motivation and discussion. — As is well known, the 
hard-sphere (HS) model is of great importance in con- 
densed matter, colloids science, and liquid state theory 
from both academic and practical points of view [l|-|3j|. 
The model has also attracted a lot of interest because it 
provides a nice example of the rare existence of nontriv- 
ial exact solutions of an integral-equation theory, namely 
the Percus- Yevick (PY) theory W for odd dimensions 

As generally expected from an approximate theory, the 
radial distribution function (RDF) provided by the PY 
integral equation suffers from thermodynamic inconsis- 
tencies; i.e., the thermodynamic quantities derived from 
the same RDF via different routes are not necessar- 
ily mutually consistent. In particular, the PY solution 
for three-dimensional one-component HSs of diameter a 
yields the following expression for the compressibility fac- 
tor Z = p/pksT (where p is the pressure, p is the number- 
density, ks is Boltzmann's constant, and T is the tem- 
perature) through the virial (or pressure) route 0-0]: 
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Here, rj = ^pc 3 is the packing fraction and the subscript 
v is used to emphasize that the result corresponds to the 
virial route. In contrast, the compressibility route yields 
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Equation is also obtained from the scaled-particle 
theory (SPT) [l6l-iT8j . The celebrated and accurate 
Carnahan-Starling (CS) [3] equation of state (EOS) is 
obtained as the simple interpolation 
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For general interaction potentials, the third conven- 
tional route to the EOS is the energy route [l|. However, 
this route is useless in the case of HSs since the internal 
energy is just that of an ideal gas, and thus it is inde- 
pendent of density. On the other hand, starting from a 
square-shoulder interaction and then taking the limit of 
vanishing shoulder width, it has been proved that the re- 
sulting HS EOS coincides exactly with the one obtained 
through the virial route, regardless of the approximation 
used [2(], [2l[ . Therefore, the energy and virial routes to 
the EOS can be considered as equivalent in the case of 
HS fluids. 



Except perhaps in the context of the SPT [161-1181] , little 
attention has been paid to a fourth route to the EOS of 
HSs: the chemical-potential route (22|- In particular, 
to the best of the author's knowledge, the possibility of 
obtaining the EOS via this route by exploiting the exact 
solution of the PY equation for HS mixtures |8| seems to 
have been overlooked. The main aim of this Letter is to 
fill this gap and derive the results 
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where /i cx is the excess chemical potential and the sub- 
script p in Eq. ([5]) denotes that the compressibility factor 
is obtained from Eq. (|4|. Equation (|5|) differs from Eqs. 
(H|) and @ in that it includes a logarithmic term and 
thus it is not purely algebraic. Nevertheless, ZpY_ fJi (r]) is 
analytic at rj = and provides well-defined values for the 
(reduced) virial coefficients b n defined by 
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Table U compares the first ten virial coefficients obtained 
from the three PY EOS, Eqs. (O, ©, and ©, with the 
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TABLE I. First ten virial coefficients b„ as obtained exactly and from several EOS related to the PY theory. 
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exact analytical (n = 2-4) and Monte Carlo (n = 5-10) 
values [23l [2D]. The interpolated coefficients obtained 
from Eq. are also included. We observe that the virial 
coefficients b 
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obtained from Eq. (J5J) are in general 
noninteger rational numbers. More explicitly, b n = 
(18-31n+15n 2 )/2n, while &l PY " u) = 2(3n-4), b { n Y ~ c) = 
(3n 2 - 3n + 2)/2, and 6^ CS) = (n + 2)(n-l). 

Interestingly enough, the virial coefficients from the 
chemical-potential route are more accurate than those 
from the virial route, although less than the ones from 
the compressibility route. This suggests the possibility 
of exploring CS-like interpolations of the form Z^irj) = 
aZpY-iiiv) + (1 — ot)ZpY-ciji) with a « 0.4. Two simple 
and convenient choices are a = | and a = j^. Thus, 

2 3 

Z^lirf) = ^PY-/i(l) + r^PY-cM 
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The values b^ c,1) = (36 - 56n + 21?i 2 + 9?i 3 )/10n and 
b { ^ 2) = (42 - 65n + 24n 2 + lln 3 )/12n obtained from 
Eqs. ([7]) and (J5J), respectively, are also displayed in Table 
HI We observe a very good agreement, even better than 
that of &,\ , with the exact values, especially in the case 



of b)l 
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In particular, b 
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is excellent. 



The superiority of Z^i and Z^ c ^ over Zcs is con- 
firmed by Fig. [TJ where the differences Zcs (v) ~ Zmd (v) j 
Z^ c ,i(v) - Zmd(v), and Z^^iv) ~ Z M b(v) (where Z M p, 
denotes molecular dynamics simulation values (25j ) are 
compared. As can be seen, both i? M c,i and ^ M c,2 deviate 
from Zmd less than Zcs over most of the stable liquid 
region. It is noteworthy that, while Z^\ predicts bet- 
ter virial coefficients than Z I1C ^ 1 the latter EOS is more 
accurate for r\ > 0.1. 

In the case of a polydisperse HS fluid characterized by 
a size distribution x(ct'), it will be proved elsewhere [26J 



that Eqs. (Q| and ([5]) are generalized to 
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where M n = da' x(a')a' n is the nth moment of the 
size distribution and in Eq. §§§ fip\(r),a) is the excess 
chemical potential of spheres of diameter a. Equation 
(1101) differs from Zpy-v and Zpy~ c for mixtures [8| only 
by the coefficient of M 2 /M|. Analogously to the one- 
component case, it is possible to construct interpolated 
EOS Z^c = aZpY-n + (1 — a)ZpY- c , with a = | or a = 
, which are more accurate than the Boubhk-Mansoori- 
Carnahan-Starling-Leland EOS (HHH. 

Once the main features and applications of the new 
PY EOS ([5]) have been discussed, the rest of the Letter 
is devoted to its derivation. 

The chemical-potential route. — Let us consider a 
d-dimensional system made of N hard spheres (i = 
1,...,N) of diameter a (the "solvent" ) plus one "solute" 
particle (i = 0) which interacts with the solvent particles 
via a HS potential of core Act with < A < 1. Thus, the 
total potential energy function is 

<S> N+1 (r N+1 ; A) = <f> N+1 (r , ri , . . . , r jv ; A) 

N 

= ^ H s ( r oj I A ) + <f>ns (nj ) , 

j=l l<i<j<N 

(11) 
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FIG. 1. (color online). Plot of Z C s(v) ~ Z M d(v) (cir- 
cles), Z^ c ,i(rj) — Zmd(ij) (triangles), and Z^ c ,2(rj) - Zmo(v) 
(squares). 



where Tij = \v,i — Yj\ and 0Hs( r ) = oo if r < er and 
otherwise. Thus, 
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where j3 = 1/ksT and 8(x) is the Heaviside step func- 
tion. We will further need the property 
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The configuration integral of the solvent+solute sys- 
tem is defined by 



Qn+i(V;X) = V-( n+1 '> Jdr Jdn... J dv N 
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where V is the volume. Note that Qjv+i(V";0) = Qn{V) 
is the configuration integral of the iV-particle system of 
solvent particles. Likewise, Qn+i{V]1) = Qn+i{V) is 
the configuration integral of a normal system of TV + 1 
identical particles. Therefore, one can write [lij 
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Qn+i(V) 



k B T ~ Q N (V) 

dX-^lnQ N+1 (V;X). (15) 

If A < j j the solute-solvent HS interaction is nonaddi- 
tive since the solute can "penetrate" the hard core of ra- 
dius \a. Provided the solvent particles are in a nonover- 
lapping configuration, the volume they exclude to the 



position of the solute particle is simply Nvd(2Xa) d , i.e., 

N 

J dr H e(r 0j — Act) = V — Nv d {2Xa) d , (16) 

where v d = (it / A) d / 2 /T(l + d/2) is the volume of a d- 
dimcnsional sphere of unit diameter. Consequently, 

Q N+1 (V; A) = [1 - (2X) d v ] Q N (V), X < i (17) 

where r\ = pvd<J d is the packing fraction of the solvent 
system, p = N/V being the number density. Equation 
([TT]) allows one to rewrite Eq. (fT5")) as 
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Qn(V) ~Q n+x {V;\) 
= - ln(l - n) - J dA A In Qiv + i(F; A). (18) 

Next, we take into account that the RDF for the solute 
particle is defined as 

y-(JV-i) 
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Application of Eq. (fT3|) into Eq. ([l"4"f then yields 
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where in the second equality g(X;r]) = g{r = Xcr + ;X,r]). 
From Eqs. (fl~7| and (|20|) we obtain the exact result 



g(A;T/)= l-(2A)V A -^' 
Let us now consider Eq. (|T8|) . Inserting Eq. ([20 
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The thermodynamic relation 



= - ln(l - 77) + d2 a i 1 I dXX^giXirj). (22) 



dp) T \dp /T 



(23) 



yields 
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Integrating both sides over density, we finally obtain 



(24) 



V) +d2 d V / 1 dAA d - 1 . 9 (A;r?) 



jnd rn pi 
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This constitutes the chemical-potential route to the EOS 
of the HS fluid. As in the virial route 



Z v ( V ) = l + 2 d - 1 V g(l;r 1 ), 



(26) 



Eq. (|25[) requires the contact value of the RDF, not the 
full spatial dependence (as required by the compressibil- 
ity route) . On the other hand, in contrast to Eq. (|26|) , Eq. 
(|25|) is "nonlocal" in the sense that it needs the knowl- 
edge of the contact value of a solute particle of diameter 
o"o = (2A — l)cr in the range < ctq < a and, moreover, 
for packing fractions rf < n. 

Equation ([25]) is formally exact. We now specialize to 
the three-dimensional case {d = 3) and consider the PY 
approximation for g(X;n), namely Q 



5py(A; rj) = 



1 



V 



2A- 1 



i- v 2(1- v y 



A 



A > 



1 



(27) 



Insertion of the above expression into the right-hand sides 
of Eqs. (J221), ([25]), and 426) finally provides Eqs. Q, ©, 
and J]}, respectively. 

The fact that Eq. ([5]) is more accurate than Eq. ([I]) can 
be explained by the following argument. Equations (|22[) . 
([24]) . and ([25]) show that the chemical-potential route is 
directly related to the integral 



Conclusion. — In summary, a hitherto hidden EOS 
for a HS fluid described by the PY liquid state theory, Eq. 
(O, has been unveiled. This new EOS from the chemical- 
potential route competes favorably with the conventional 
one from the virial route by the reasons outlined above. 
Thus, at least in the framework of the PY theory, the 
chemical-potential route should be placed on the same 
footing as the standard virial and compressibility routes. 
Apart from its intrinsic academic and pedagogical inter- 
est, the new EOS has a practical impact. For instance, 
the chemical-potential and compressibility routes allow 
for the construction of interpolation proposals, Eqs. ([7]) 
and ([5]), which are more accurate than the widely used CS 
EOS. Moreover, the use of Eq. ^ for HS fluids and Eq. 

for more general systems can be very helpful for the 
construction of accurate EOS. Extensions of this work to 
sticky hard spheres and to hyperspheres d, Gil 13 1 
are planned. 
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dAA g(A;jj), 



(28) 



which is proportional to a (weighted) average of the con- 
tact value g(X; n) in the range \ < A < 1. Since both the 
contact value g(X;r]) and its first derivative dg(X;rj)/dX 
at A = | are given exactly by the PY equation [compare 
Eqs. (]21l) and ([27]) ]. it seems reasonable that the "aver- 
age" value ([28]) is better estimated than the end point at 
A = 1 by the PY approximation. 

Thus far, all the results have been specialized to HS 
systems. In the more general case of particles interacting 
through a potential <f>{r) , one can still single out a particle 
i = which interacts with the rest via a potential (j)(r; A) 
such that <p(r; 0) — and (f>(r; 1) = <t>(r). Proceeding in 
a similar way as before, one arrives at [261 ] 
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In this equation the A-protocol (j>(r; X) remains arbi- 
trary. If 4>{r) is not a singular potential, an obvious 
choice is (j)(r; A) = X<f)(r) [29j. However, this possibil- 
ity is ill-defined if, as happens with ^Hs( r ), the poten- 
tial diverges over a finite range. In that case, an ad- 
equate choice is </>(r; A) = (f>(r/X). Using the identity 



A 0(r;A) = -k B Te^ r ^dxe-^ r ^ 



particularizing to 
yields Eq. ([2"2"]> . 
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